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ABSTRACT 

GRB 090426 is a short duration burst detected by Swift (Tgo ^ 1.28 s in the observer 
frame, and Tgo ~ 0.33 s in the burst frame at z = 2.609). Its host galaxy properties and 
some 7-ray related correlations arc analogous to those seen in long duration GRBs, 
which are believed to be of a massive-star origin (so-called Type II GRBs) . Wc present 
the results of its early optical observations with the 0.8-m TNT telescope at Xinglong 
observatory, and the 1-m LOAO telescope at Mt. Lemmon Optical Astronomy Ob- 
servatory in Arizona. Our well-sampled optical afterglow lightcurve covers from ~ 90 
seconds to ~ 10"^ seconds post the GRB trigger. It shows two shallow decay episodes 
that are likely due to energy injection, which end at ~ 230 seconds and ^ 7100 seconds, 
respectively. The decay slopes post the injection phases are consistent with each other 
(a ~ 1.22). The X-ray afterglow lightcurve appears to trace the optical, although the 
second energy injection phase was missed due to visibility constraints introduced by 
the Swift orbit. The X-ray spectral index is Px ^ 1-0 without temporal evolution. Its 
decay slope is consistent with the prediction of the forward shock model. Both X-ray 
and optical emission is consistent with being in the same spectral regime above the 
cooling frequency {vc)- The fact that Vc is below the optical band from the very early 
epoch of the observation provides a constraint on the burst environment, which is 
similar to that seen in classical long duration GRBs. We therefore suggest that death 
of a massive star is the possible progenitor of this short burst. 

Key words: gamma rays: bursts (individual: GRB 090426) — gamma rays: observa- 
tions. 



1 INTRODUCTION 

Cosmic gamma-ray bursts (GRBs) are classified into two 
classes with a separation at the observed burst duration 
of Tgo ~ 2 seconds based on CGRO/BATSE observations 
(Kouveliotou et al. 1993). The afterglow and host galaxy 
properties of long GRBs, especially the detections of several 
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long GRB-supernova associations (e.g. Galama et al. 1998; 
Hjorth et al. 2003; Soderberg et al. 2004; Campana et al. 
2006), suggest that they are mostly likely related to deaths 
of massive stars. The "coUapsar" model has been widely rec- 
ognized as the standard scenario for long GRBs (Woosley 
1993; Paczyhski 1998; Zhang & Meszaros 2004; Piran 2004; 
Woosley & Bloom 2006). In the Swift era, the afterglows 
and host galaxies of some short GRBs have been identified 
(Gehrels et al. 2005; Villasenor et al. 2005; Fox et al. 2005; 
Berger et al. 2005; McGlynn et al. 2008). Some short GRBs 
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are found to be associated with nearby early-type galaxies 
with little star formation. Some others are located in late- 
type star forming galaxies (e.g. Fox et al. 2005; Fong et al. 
2010), some of which are at high redshifts (Levan et al. 2006; 
Berger et al. 2007) . No short-duration GRB was found to be 
associated with a supernova (Kann et al. 2008; Zhang et al. 
2009, and references therein). All these seem to favor the 
idea that short GRBs are from mergers of two compact stel- 
lar objects (Eichler et al. 1989; Narayan et al. 1992). 

Later Swift observations have revealed that the 
short /long separation is not sufficiently clean to differentiate 
their physical origins. This led to introduction of the physi- 
cal classification scheme Type If (collapses of massive stars) 
vs. Type I (putatively identified as mergers) GRBs (Zhang 
et al. 2007a; Kann et al. 2007, 2008; Zhang et al. 2009; Lv et 
al. 2010). Some convincing Type I GRBs have long, soft "ex- 
tended emission" (Barthelmy et al. 2005; Norris et al. 2006; 
Lin et al. 2008; Zhang et al. 2009; Perley et al. 2009), mak- 
ing their Tgo "long". The non-detection of any supernovae 
associated with the nearby long GRBs 060614 and 060505 
(Gehrels et al. 2006; Gal- Yam et al. 2006; Fynbo et al. 2006) 
casted doubts on the Type II origin for these long GRBs. 
Some long GRBs have rest-frame durations shorter than 2 
seconds (Levan et al. 2007). Observations of two intrinsically 
short-duration, high-z GRBs 080913 {z = 6.7; Greiner et al. 
2009a) and 090423 {z = 8.3; Tanvir et al. 2009; Salvaterra et 
al. 2009) suggest that they share a lot of common properties 
with long GRBs, and most likely have a massive star progen- 
itor (Zhang et al. 2009; Lin et al. 2009; Levesque et al. 2010a; 
Belczynski et al. 2010). Zhang et al. (2009) and Virgih et al. 
(2009) suggested that some (or even most) short duration 
GRBs are probably not produced via compact star mergers 
(Type I), but are likely related to massive stars (Type II). 
Based on the observed gamma-ray energy and peak energy 
of the i^/y spectrum of prompt gamma-ray emission, Lv et 
al. (2010) defined a parameter e = Eiso/Ep l, and proposed 
a new empirical classification scheme that is found to bet- 
ter match the physically-motivated Type II/I classification 
scheme. They showed that the typical Type II GRBs are in 
the high-e group, in contrast to the typical type I GRBs, 
which belong to the low-e group. 

Another striking case that poses a challenge to the con- 
ventional long vs. short GRB classification scheme is GRB 
090426. This burst has an observed Tgo = 1.29±0.09 seconds 
in the Swift BAT band (Sato et al. 2009), corresponding to 
a burst rest-frame duration of 0.33 s at redshift z — 2.609 
(Levesque et al. 2010a). Phenomenologically, in terms of du- 
ration alone, it is unambiguously within the range of clas- 
sical short-type GRBs in both the observer frame and the 
burst rest frame. On the other hand, both the host galaxy 
properties (Levesque et al. 2010a) and the spectral energy 
properties (Antonelli et al. 2009) suggest that it is mostly 
likely a Type II GRB, i.e., related to core collapse of a mas- 
sive star. With the new classification method proposed by Lv 
et al. (2010), this event is also well classified into the high-e 
group, which is where all known Type II GRBs belongs to. 

Multi-wavelength afterglows are essential for revealing 
the burst environment, and hence can serve as a probe of the 
GRB progenitor. In this paper, we report our observations of 
the early optical afterglow for GRB 090426 using the TNT 
telescope at Xinglong Observatory and the LOAO telescope 
at Mt. Lemmon Optical Astronomy Observatory in Arizona. 



We use the optical and X-ray afterglow data to explore the 
nature of this event. Our observations are reported in Sec- 
tion 2. A joint optical and X-ray data analysis is present in 
Section 3. Conclusions and discussion are present in Section 
4. The notation f^, oc t -"u-'' is used throughout the paper, 
where is the spectral flux density at the frequency ly. 



2 OBSERVATIONS 

GRB 090416 was detected by Swift Burst Alert Telescope 
(BAT) at 12:48:47 UT on 2009 April 26 (Cummings et al. 
2009). Its duration is Tgo = 1-28 ± 0.09 sec in 15-350 KeV. 
The Swift X-Ray Telescope (XRT) began to observe the 
burst since 84.6 sec after the GRB trigger. At 89 sec after 
the trigger. The Swift UVOT began to observe the burst and 
reported an optical counterpart with a brightness of about 
17.5 mag in the white band. The optical afterglow was con- 
firmed by the Xinglong TNT telescope (Xin et al. 2009) and 
other follow-up observations (e.g. Im et al. 2009). A redshift 
of 2.609 was determined by Levesque et al. (2009c) using 
the Keck telescope, which was confirmed by the ESO VLT 
observation (Thoene et al. 2009). The time-integrated 7-ray 
spectrum is well fit by a power-law with a photon index F 
of 1.93 (Sato et al. 2009) (/J = F - 1), which roughly corre- 
sponds to an estimated spectral peak energy (in the observer 
frame) Ep^obs ~ 45 keV using an empirical relation between 
Ep and the power-law photon index of the BAT spectrum 
(Zhang et al. 2007a; Sakamoto et al. 2009). 



2.1 Optical Observation and Data Reduction 

We carried out a follow-up observation campaign of GRB 
090426 using the TNT (0.8-m Tsinghua University - Na- 
tional Astronomical Observatory of China Telescope) at 
Xinglong Observatory, under the framework of East-Asia 
GRB Follow-up Observation Network (EAFON, Urata et al. 
2003; 2005). TNT is equipped with a PI 1300 x 1340 CCD 
and filters in the standard Johnson Bessel system. Its field of 
view is 11.4 x 11.4 arcmin, yielding a 0.5 arcsec pixel scale. 
A custom-designed automation system has been developed 
for the GRB follow-up observations (Zheng et al. 2008). 

The observation of the optical transient (OT) of GRB 
090426 was carried out with TNT at 86 seconds post the 
Swift/BAT trigger, which is slightly earlier than the begin- 
ning observation of XRT and UVOT on-board Swift. A new 
fading source was discovered and confirmed as the OT of 
the burst (Xin et al. 2009). The coordinates of the OT are 
consistent with that of UVOT (Cummings et al. 2009). The 
W (white), R and V-band images were obtained in 86—519 
seconds, 570 — 1513 seconds, and 1907 — 10748 seconds post 
the GRB trigger, respectively. 

The 1-m telescope (Han et al. 2005) is located at Mt. 
Lemmon Optical Astronomy Observatory LOAO in Arizona 
operated by the Korea Astronomy Space Science Institute. 
The observation of GRB 090426 was carried out at about 
16.3 hours after the burst {t ~ 60 ks) in the R filter. The 
OT was not detected and only upper limits were obtained. 

Data reduction was carried out following the stan- 
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dard routine in package, including bias and flat- 

field corrections. Dark correction was not performed since 
the temperature of our CCD was cooled down to — 110°C. 
Point spread function (PSF) photometry was applied via 
the DAOPHOT task in the IRAF package to obtain the in- 
strumental magnitudes. During the reduction, some frames 
were combined in order to increase the signal-to-noise ratio 
(S/N). In the calibration and analysis, the white band was 
treated as the R band (Xin et al. 2010). Absolute calibration 
was performed using the Sloan Digital Sky Survey (SDSS, 
Adelman- McCarthy et al. 2008), with conversion of SDSS 
to Johnson- Cousins systenQ- The data of GRB 090426 ob- 
tained by TNT and LOAO are reported in Table. [T] 

2.2 Swift/XRT X-ray Afterglow Data Reduction 

The SwiftfKKT lightcurve and spectrum are extracted from 
the UK Swift Science Data Centre at the University of Le- 
icester (Evans et al. 20090- We fit the X-ray spectrum with 
the Xspec package. The time-integrated X-ray spectrum is 
well fit by an absorbed power-law model, with a photon 
power-law index F = 2.00 ± 0.06 from the PC mode data. 
No significant host Nh excess over the Galactic value is de- 
tected. 



3 OPTICAL AND X-RAY AFTERGLOW 
JOINT ANALYSIS 

3.1 Temporal Analysis 

With the X-ray spectral index (/3x ~ 1-00 ± 0.06), we first 
derive the 1 keV lightcurve from the XRT data. It is shown 
in Fig. [T] Next, we convert the extinction-corrected magni- 
tudes of the optical afterglow into energy fluxes. Levesque 
et al. (2010a) reported Av ~ 0.4 for the GRB host galaxy, 
and the extinction by the MilkyWay Galaxy in the burst di- 
rection is Eb-v = 0.017 (Schlegel et al.l998). For the host 
galaxy, the transformation from Ay to Ar is nearly inde- 
pendent of any known type of extinction laws (MW, LMC, 
SMC), and we take Ar ~ 0.32 for the host galaxy. The Milky 
Way extinction corresponds to Ar = 0.046 and Ay = 0.057. 
After extinction corrections, we find that the uF,^ fluxes in 
both the optical and the X-ray bands are almost the same. 
This indicates a flat vF^ spectrum from the optical to the 
X-ray bands, i.e., Pox ~ 1.0. This is consistent with the 
observed X-ray spectral index. This flat uF,j spectrum also 
means that we do not need to calibrate the observed opti- 
cal fluxes to a given band. The extinction-corrected optical 
energy flux lightcurve is also shown in Fig. [1] 

As shown in Fig. [T] the afterglow lightcurves in the op- 
tical and X-ray bands trace with each other, suggesting that 
not only they are of the similar (external shock) origin, but 
they also belong to the same spectral regime. At early epochs 
{t < 2200 seconds post the GRB trigger), the optical and X- 
ray afterglow lightcurves show a smooth shallow-to-normal 

1 IRAF is distributed by NOAO, wliicli is operated by AURA, 
Inc., under cooperative agreement with NSF. 

^ http:/ /www. sdss.org/dr6/algorithms/sdssUBVRITransform. html 
#Lupton2005 

http://www.swift.ac.uk/rcsults.shtml 



Table 1. Optical Afterglow Photometry Log of GRB 090426. The 
reference time Tq is Swift BAT burst trigger time. All Data are 
not corrected for the Galactic extinction (which is Er-v = 0.017, 
Schlegel et al.l998). "Merr" means the uncertainty of magnitude. 
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Table 2. Optical lightcurve fits with a smooth broken power-law 
model for two cpoches. 



Interval Fq *p(s) ck2 



86-2200S 3.74±0.32 227±27 0.26±0.07 1.22±0.04 
> 2200 s 1.77±0.67 7123±3620 0.19±0.41 1.22 



In units of xlO erg cm ^ s 
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Figure 1. The extinction-corrected optical light curve of GRB 
090426 (sohd dots) and the best fit with a smooth broken power- 
law model (solid line) for two epoches before and after 2200 
seconds post the GRB trigger. The i?-band data observed by 
Rumyantsev et al. (2009) is marked with squares. The XRT 
lightcurve (open dots) is also shown for comparison. 



transition in the decaying behavior, with small Bickerings in 
the X-ray band. In the time interval 2200 — 5500 seconds, the 
brightness of the OT was almost constant, fading again after 
t > 5500 seconds. The OT was also detected by Rumyantsev 
et al. (2009) at 0.2342 days (20200 seconds) post the GRB 
trigger. It had faded to i? = 21.3 ± 0.2, indicating a decay 
slope of ~ 1.2 in the time interval oit> 5500 seconds. Due 
to the orbital constraints of the Swift satellite, there is no 
XRT observation in the time interval [2200, 5500] seconds 
when the optical lightcurve features a plateau. The fading 
behavior of X-rays in the time interval [5500, 7700] seconds 
is similar to that of the optical afterglow. At t ~ 0.2342 day, 
the detected X-ray behavior is also consistent with that of 
the optical emission. These results suggest that the tempo- 
ral behavior of both optical and X-ray afterglows could be 
the same. We thus perform a temporal analysis on the well- 
sampled optical lightcurve only. The optical lightcurves in 
the time intervals both before and after 2200 seconds can be 
well fit with a smooth broken power-law. 



F = Fo 



+ 



(1) 



Our best-fit parameters are summarized in Table. (2] We find 
that the decay slopes post the two breaks are similar, with 
a value ~ 1.22. The late time X-ray data at ~ 4 x 10^ s is 
also consistent with the predicted behavior of such a decay 
slope in the X-ray band. 



3.2 Data Confronted with the Forward Shock 
Models 

As shown above, the temporal behaviors of the optical and 
X-ray data are consistent with being achromatic, and the 
decay slopes after the two breaks observed in the optical 
lightcurve are well consistent with the prediction of the for- 
ward shock models. The fact that Oo = is consistent 



with the forward shock model in the same spectra regime 
(see also Urata et al. (2007) for a more general discussion of 
the Oo — Qa; relation for the forward shock models). We find 
no spectral evolution across the two breaks from the X-ray 
data. The derived I5x are 1.09 ± 0.15 and 1.03 ± 0.10, re- 
spectively, for the integrated X-ray spectra a±t < 2200 s and 
t > 2200 s. Inspecting the spectral index Px and the tempo- 
ral decay index ax in the normal decay segment, one finds 
good agreement between data and the forward shock model 
closure relation in the spectral regime v > max^Vm, Vc), i.e. 
a = {SPx — 1)/2 = 1.14±0.23, where i/m and Vc are the typi- 
cal and cooling frequencies of synchrotron radiation, respec- 
tively. Since both optical and X-ray emissions are consistent 
with the forward shock origin in the normal decay phase, 
one can naturally attribute the two shallow decay segments 
to two epochs of energy injection into the blastwave (e.g., 
Dai & Lu 1998; Zhang & Meszaros 2001; Zhang et al. 2006; 
Liang et al. 2007). 



3.3 Constraints on the Burst Environment 

The circumburst environment is critical to understand the 
nature of a GRB. In the literature, usually two types of 
medium are discussed, namely, a constant density medium 
relevant for interstellar medium (ISM) or a stratified stellar 
wind with a density profile n oc r~^. For the spectral regime 
identified for GRB 090426, i.e. v > max(z^„i, Vc), unfor- 
tunately the observed flux does not depend on the medium 
density. Consequently, one can not distinguish the two types 
of medium. However, since a wind model would undoubt- 
edly point towards a massive star progenitor, we only focus 
on the constant density case. According to the analysis in 
section 3.2, Vc should be below the optical band at very 
early epochs. This would give an interesting constraint on 
the medium density. 

In the constant density case, the typical synchrotron 
emission frequency, the cooling frequency and the peak spec- 
tral flux density are (Sari et al. 1998; coefficients taken from 
Yost et al. 2003; Zhang et al. 2007b): 



= 3.3x10'^ Hz 



p-2 
p-1 



1/2,1/2 2 



pl/2 ,-3/2 

(2) 



6.3x10^^ Y{z(l + z)-^'\l + Y)-\-^%E-^]l^n~h-^/^ 

(3) 



1.6 mjy(l -I- z)Z)28^'^B -2-'^^>52n 



,1/2 



(4) 



where D is the luminosity distance in units of 10 cm, /p 
is a function of p (/p ~ 1 for p = 2, Zhang et al. 2007b), 
and td is the observer's time in unit of days. The convention 
Qn = Q/IO" is adopted in cgs units. The Inverse Compton 
scattering parameter 



y = [-1 + (IH- 477i772eeAs)'/']/2 



(5) 



where 7?i = min[l, (i^c/i'm)'^"*''^^] ~ 1 for p ~ 2.0 (Sari & 
Esin 2001), and 772 ^ 1 is a correction factor introduced by 
the Klein-Nishina effect. We take 772 = 1 in our analysis, so 
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that the _R-band energy flux reads 



4 DISCUSSION AND CONCLUSIONS 



where 



F 1/2 (p-l)/2 (2-p)/2 



= 2.67 X 10""/p ergs s"' cm-^D^sil + 2^"+^^^^ 
X (1 + *<-\<S'/^tr'''''''-«''-''''(%) 



7.45 X 10" 



p-1 



{p-l)/2 



(7) 



Since p = 2/3 = 2.18 ± 0.30, we have fp = 6.04 x 10"^. One 
can then derive Ek using the data at any time td (Zhang at 
al. 2007b): 



E 



K,52 = 



1 4/(p+2) 



1.61 X 10"^-^ ergs s~^ cm"^ 



l,(3p-2)/(p+2) 



X (i + y)^/(*'+2)4^72^/(''+^) 

V 4{l-p)/{p+2) 2(p-2)/(p+2) 



(8) 



The analysis above is valid for p > 2. From the X-ray data, 
we find the p is slightly larger than 2. In order to make a 
rough estimate, we take p ~ 2 to simplify Eq. [S] as 



i^rF^ 



1.61 X 10~^3 ergs s"i cm"^ 1 + z 
The cooling frequency Uc is given by 

i^rFuj^ 



6.3 X 10'" Hz 



6r/_i(i+F)z?i,28. 

(9) 



-1/2 



1.61 X 10^^^ ergs s^^ cm^^ 



(10) 



At t ~ 100 seconds, one has vrF^j^ = 4.42 x 10 erg cm ^ 
, so that 



Vc = 1.57 X 10^^ Hz(l + y)" 



'=S,-2'=e,-l' 



(11) 



The requirement of ur > Vc eA, t ^ 100 seconds then gives 



n > 354s(ee, es) 



(12) 



where s(ee, es) = £^^^2^^^1(1-1-^)"^''^. This already points 
towards a large medium density. Noticing that s would de- 
crease when es increases, we consider an extreme case of en- 
ergy equipartition among radiating electrons, magnetic field 
and baryons, i.e., — 1/3, es = 1/3. In this case, we still 
get n > 11.2 cm~'^. Such a dense medium is inconsistent 
with the naive expectation of a compact star merger pro- 
genitor, which tend to occur in low-density medium as the 
compact binary escapes from the star forming region due to 
the natal kicks during the births of the two neutron stars. 
The high density, on the other hand, is consistent with that 
expected for a massive star that was born in a high density 
star forming region such as a molecular cloud (e.g. Chevalier 
et al. 2004). 

Recall that a stellar wind medium also suggests a mas- 
sive star connection, we therefore conclude that the high 
ISM density would favor a burst born in a region of ongoing 
star formation, and hence favor a massive star progenitor 
(or a Type II GRB) for GRB 090426. 



We have presented the early optical afterglow observations of 
GRB 090426 with TNT and LOAO. Our well-sampled opti- 
cal afterglow lightcurve from ~ 90 seconds to ~ 10* seconds 
post the GRB trigger seems to exhibit two energy injection 
phases that ended at ~ 230 seconds and ~ 7100 seconds, 
respectively. The temporal behavior of the observed X-rays 
is consistent with the optical emission. We show that both 
the optical and X-ray afterglows are consistent with the for- 
ward shock model in the spectral regime u > max{um, Vc)- 
Although the medium type (ISM vs. wind) cannot be dis- 
tinguished, we can make a case of the massive star (Type 
II) origin of the GRB. For a constant density medium, the 
required medium density is found very high (Eq. [12]), in- 
consistent with the expectation of a compact star merger 
(Type I) progenitor. A massive star progenitor forming in 
a high-density star forming region (e.g. molecular cloud) is 
consistent with the data. If the circumburst medium is a 
stellar wind, then it also directly points towards a massive 
star progenitor. We therefore conclude that the afterglow 
data of the short GRB 09046 strongly suggests that it is a 
Type II GRB. 

GRB 090426 is of great interest because of its short 
duration (Tgo ~ 0.33 seconds in the burst frame) and high 
redshift [z = 2.609). Its redshift significantly exceeds the 
previous spectroscopically confirmed high redshifts of short 
GRBs, e.g. GRB 070714B {z = 0.904, Graham et al. 2009), 
GRB 051121A (2 = 0.546, Soderberg et al. 2006b), GRB 
070429(2 = 0.902, Genko et al. 2008) and GRB 090510 
{z — 0.903, Rau et al. 2009). Its host galaxy is a blue, lumi- 
nous, and star-forming galaxy (Levesque et al. 2010a), sim- 
ilar to the host galaxies of typical Type II GRBs (Levesque 
et al. 2010b). The spectral energy correlation is also consis- 
tent with being a Type II GRB. Here we present a third 
argument in favor of its Type II origin, based on early af- 
terglow observation and modeling. It is interesting to note 
that with the new classification method proposed by Lv 
et al. (2010), GRB 090426 is well grouped into Type II 
GRBs, in contrast to other short GRBs such as GRB 051221 
and GRB 070714B. We also noted that GRB 090426 shares 
some similar properties with two other Type II bursts, GRB 
040924 and GRB 050416A. The first was a short-duration 
(1.2-1.5 sec), soft spectrum (Huang et al. 2005) and a SN 
1998bw-like supernova-associated burst (Soderberg et al. 
2006a; Wiersema et al. 2008). The second burst also ex- 
hibited a soft spectrum and a short duration in the rest 
frame (Tgo = 2.4s, z = 0.6528), which is associated with 
a SN 1998bw-like supernova (Soderberg et al. 2007), and is 
located in a circumburst medium with a large density vari- 
ation. 

Finally, GRB 090426 has reinforced the conclusion that 
burst duration alone cannot be used to judge the physical 
nature of a GRB, as discussed by Donaghy et al. (2006) and 
Zhang et al. (2007aj3; Levan et al. (20070, and Zhang et 



* Both teams proposed that some long GRBs can be physically 
associated with compact star mergers, but did not discuss that 
short bursts can be physical related to massive star core collapses. 
^ These authors pointed out that some long GRBs are rest-frame 
short. 
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Multi-wavelength observational campaigns are 
essential to unveil the physical nature of GRBs. 
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